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I. INTRODUCTION 

Although the "two-body problem" cannot be solved 
analytically in the general theory of relativity, at least 
two approximation methods can be used to tackle it. 
The first one dates back to Einstein's 1915 calculation 
of the relativistic perihelion advance of Mercury's orbit 
[1] , and is based on a perturbative treatment in powers of 
the ratio v/c between the binary's relative velocity v (in 
the center-of-mass frame) and the vacuum speed of light 
c. At the lowest order of approximation, this approach 
gives back the well-known Newtonian solution, and is 
therefore dubbed "post-Newtonian" (PN) expansion; see 
e.g. Ref. [2] and references therein. Currently, the two- 
body dynamics of non-spinning compact objects is known 
through 3PN order 1 [3, 4], and the gravitational- wave 
fluxes of energy and angular momentum through 3.5PN 
[5-7] and 3PN [8] orders for circular and eccentric or- 
bits, respectively. Spin effects have also been computed 
in the dynamics and gravitational radiation [9-11] using 
both Hamiltonian [12-16] and Lagrangian [9, 17] formu- 
lations. The effective-field-theory approach applied to 
gravity [18] has confirmed some of these results [19, 20], 
and has pushed the calculations to higher PN orders for 
spinning bodies [20-27]. 

The second approximation method also dates back to 
1915, namely to Karl Schwarzschild's famous wartime 
calculation of the gravitational field of a spherically sym- 
metric body [28, 29]. While non-rotating black holes are 



As usual we refer to nPN as the order corresponding to terms 
0(c~ 2n ) with respect to the Newtonian acceleration in the equa- 
tions of motion, or with respect to the quadrupole formula in the 
radiation field. 



described by the Schwarzschild metric, rotating black 
holes are represented by the Kerr solution [30]. The 
motion of test masses in the Schwarzschild or Kerr ge- 
ometries is naturally described by the geodesic equa- 
tions, which are valid for arbitrarily high values of the 
orbital velocity v (i.e. the geodesic equations formally 
include all PN corrections). If the finite mass of the 
particle and its backreaction on the background geom- 
etry are taken into account, the orbits will deviate from 
geodesic motion under the effect of the gravitational self- 
force (GSF) [31]. More formally, Schwarzschild or Kerr 
geodesies can be seen as the orbital motion of the binary 
at the zeroth order of approximation in the mass ratio, 
while at first order the two-body dynamics is regulated 
by the so-called "MiSaTaQuWa" equation [32, 33]. The 
GSF can be split into a dissipative component related to 
gravitational-wave emission — which is described by the 
Regge- Wheeler and Zerilli equations in a Schwarzschild 
background [34, 35], and by the Teukolsky equation in 
Kerr [36, 37] — , and a conservative component respon- 
sible for secular effects such as the periastron advance 
[38-40]. 

A notable event in the history of the general relativis- 
tic two-body problem took place almost a century after 
Einstein's and Schwarzschild's early work, in 2005, with 
the first numerical simulations of the inspiral, merger, 
and ringdown of a system of two non-spinning black holes 
[41-43]. While these results constitute a truly remarkable 
achievement, current "state-of-the-art" numerical rela- 
tivity (NR) simulations are still too time-consuming to 
provide gravitational waveforms covering the whole pa- 
rameter space of binary black-hole systems, especially for 
small mass ratios [44] and black holes with large spins 
[45,46]. 

A semi-analytical approach that is flexible enough to 
incorporate information from both PN expansions and 
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black-hole perturbations, as well as from NR simulations, 
is the effective-one-body (EOB) method [47]. The basic 
idea behind this construction is to map the orbital dy- 
namics of an arbitrary mass-ratio compact binary sys- 
tem onto that of a test particle in a suitable background 
spacetime. In order for the EOB model to have the 
correct test-particle limit, this effective background met- 
ric, g 1 ^, must clearly reduce to that of a Schwarzschild 
black hole when one of the masses goes to zero (for non- 
spinning binaries). In addition, such a mapping is known 
to exist for any mass ratio at the Newtonian level, be- 
cause in Newtonian gravity one can always map the dy- 
namics of a binary system with masses mi and mi onto 
the motion of an effective particle with reduced mass 
fi = m\mij (m\ + m-i) around a body with total mass 
M = mi + mi. It is therefore natural to try to achieve 
the EOB mapping by considering an effective particle 
with mass /i moving in a time-independent and spheri- 
cally symmetric deformed Schwarzschild spacetime with 
total mass M, 2 

ds 2 cS = gtf (r; v) di 2 + gf r {r; v) dr 2 + r 2 df> 2 , (1.1) 

where the deformation is regulated by the symmetric 
mass ratio v = n/M. In the test-particle limit v — > 0, 
we recover (by construction) the Schwarzschild result 
5 t c f = -l/g°r = -1 + 2M/r. Beyond that limit, the 
^-dependence of the EOB potentials g*j.f and g°r encodes 
crucial information about the dynamics of the real binary 
system. 

Buonanno and Damour [47] showed that, for spinless 
binaries, the EOB mapping can be achieved not only at 
the Newtonian level, but also at the 1PN and 2PN orders, 
obtaining the following relation between the Hamiltonian 
-ffrcai of the real binary system and the Hamiltonian H c g 
of the effective particle: 



H Ieal =M Jl + 2v(^-lj. (1.2) 

Remarkably this formula coincides with that introduced 
by Brezin, Itzykson and Zinn- Justin [48] in quantum elec- 
trodynamics to map the one-body relativistic Balmer for- 
mula with the two-body one; for example it can relate 
some of the energy levels of positronium (an equal-mass 
system comprised of an electron and an anti-electron, de- 
scribed by the real two-body Hamiltonian Hreai) to those 
of the hydrogen atom (an "extreme" mass-ratio system 
described by the effective Hamiltonian H c g). 

Since the original paper [47], the EOB mapping has 
been extended to 3PN order for non-spinning bina- 
ries [49], and shown to exist also for spinning binaries, 
through 3.5PN order in the spin-orbit terms, and 2PN or- 
der in the spin-spin terms [50-55] . Furthermore, the EOB 
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construction has grown to include a model for the grav- 
itational waveforms [56-60], allowing detailed compar- 
isons (and calibrations of the EOB model's unknown pa- 
rameters) with NR waveforms for non-spinning and spin- 
ning comparable-mass binaries [61-71], as well as with 
Regge-Wheeler-Zerilli [57, 72-75] and Teukolsky wave- 
forms [60, 76-78] for small and extreme mass-ratios. 

More recently, information coming from GSF calcula- 
tions has started to be included into the EOB model. 
References [53, 71, 79] used the frequency shift of the 
Schwarzschild innermost stable circular orbit (ISCO) in- 
duced by the conservative GSF, as recently calculated 
by Barack and Sago [39] (see also Ref. [80]), to constrain 
some unknown parameters entering the g^f component of 
the EOB effective metric (1.1) for spinless binaries, and 
regulating the unknown higher PN orders. Ref. [79] also 
suggested using GSF data to determine a certain combi- 
nation pgp of the free functions parametrizing the g^f and 
g°r components of the EOB effective metric, at linear or- 
der in the mass ratio. That suggestion was then realized 
in Ref. [81], which calculated the strong- held behavior of 
Psf using the GSF contribution to the periastron advance 
for quasi-circular orbits in a Schwarzschild background. 

Besides the ISCO frequency shift and the periastron 
advance, a third physically meaningful (i.e. coordinate 
invariant) perturbative result that has been obtained 
within the GSF framework is the conservative effect of 
the self-force on the "redshift observable" z = —u a K a , 
u a being the particle's four-velocity and K a the helical 
Killing vector of the perturbed Schwarzschild geometry 
[38, 82-85]. Refs. [38, 84, 85] (see also [86]) demonstrated 
a very good agreement between the numerical GSF result 
and the PN prediction. 

The quantity z measures the redshift of light rays emit- 
ted from the particle, and received far away from the 
binary system, along the helical symmetry axis perpen- 
dicular to the orbital plane [38]. In the companion pa- 
per [87] , we put forward a different interpretation for the 
redshift observable. Building on the hrst law of binary 
point-particle mechanics recently established in Ref. [88] , 
we relate z to the binding energy of a binary system 
of non-spinning compact objects through next-to-leading 
order in the mass ratio. Using the numerical results of 
Refs. [38, 82-85], in which the GSF contribution to the 
redshift observable of a non-spinning particle moving on 
a circular orbit around a Schwarzschild black hole was 
calculated, we derive an explicit expression for the bind- 
ing energy. 

In this paper, we use this new result for the binding 
energy, assuming that the mapping (1.2) holds at all PN 
orders, to derive the exact expression of the g^f compo- 
nent of the EOB effective metric (1.1) for non-spinning 
binaries, through linear order in the mass ratio. The re- 
sulting expression exactly reproduces the ISCO frequency 
shift calculated by Barack and Sago, and goes far beyond 
the results of Refs. [53, 71, 79] that could only constrain 
certain combinations of unknown parameters appearing 
in the EOB potential g^f. 
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Moreover, thanks to the exact knowledge (at least 
through linear order in the mass ratio) of g^f, we are able 
to use the GSF results of Ref. [79, 81] for the periastron 
advance to determine the EOB metric component g^r 
exactly (through linear order in the mass ratio) for non- 
spinning binaries. This constitutes significant advance 
over earlier results that could only constrain combina- 
tions of unknown functions entering g^f and gf£. Our 
results therefore completely determine the EOB effective 
metric (1.1) for a system of non-spinning compact ob- 
jects, at first order in the mass ratio. 

This paper is organized as follows. In Sec. II, after 
briefly reviewing the EOB effective metric and Hamil- 
tonian dynamics, we use the binding energy computed 
in Ref. [87] to derive the exact correction to gff that is 
linear in the mass ratio. Moreover, employing the re- 
sults of the periastron advance for circular orbits from 
Refs. [79, 81], we also derive the exact term linear in 
the mass ratio in g^. In Sec. Ill we use the binding en- 
ergy computed through 6PN order in Ref. [88], together 
with several constraints among the EOB metric coeffi- 
cients derived in Refs. [79, 81], to compute <7|f and g°r 
through 6PN and 5PN orders, respectively, at linear or- 
der in the mass ratio. Finally, in Sec. IV we discuss the 
main results of the paper and comment on future work. 
The structure of the EOB Hamiltonian used in the rest 
of this paper is detailed in an Appendix. 



II. SELF-FORCE CONTRIBUTIONS TO THE 
EOB POTENTIALS 

Within the EOB framework, the real Hamiltonian 
-ff r cai encoding the orbital dynamics of two non-spinning 
compact objects is mapped to an effective Hamiltonian 
H e s describing a test particle of mass /i = m\mil {m\ + 
m<i) moving in a deformed Schwarzschild metric of mass 
M = mi + ?7i2- The deformation is regulated by the bi- 
nary's symmetric mass ratio v = fi/M, and disappears in 
the test-particle limit v — > 0. The EOB effective metric 
reads [47] 

ds 2 cS = -A(r) At 2 + B(r) dr 2 + r 2 dil 2 , (2.1) 

where the potentials A = —g^f and B = g™p are known 
through 3PN order [47, 49]. We find it more convenient 
to work with the potential D = (AB) -1 , so that 

A{u) = 1 - 2u + 2vu 3 + f y - ^tt 2 ^ vu 4 + 0{u 5 ) , 

(2.2a) 

D(u) = 1 + Qv u 2 + (52 - &v) v u 3 + 0(u 4 ) , (2.2b) 

where u = M/r denotes the inverse Schwarzschild-like 
EOB radial coordinate. In the test-particle limit v — » 0, 
we recover (by construction) the Schwarzschild results 
A(u) =l-2u and D(u) = 1. Through 3PN order, the 



effective Hamiltonian is given by [47, 49] 

H 2 S (r,p) = p 2 A(r) [l + p 2 + (B(r)- 1 - l) (n ■ p) 2 

+ 2{A - Zv) v u 2 {n-pf} , (2.3) 

where we introduced the reduced conjugate momentum 
p = p/fi and the unit vector n — r/r. Finally, the so- 
called improved real (or EOB) Hamiltonian reads [47] 

Heob = HlZ Iovcd = M ^l + 2^^£-l) . (2.4) 

When extending the EOB Hamiltonian (2.4) to higher 
PN orders, one needs to modify the effective dynamics 
with terms depending on higher-order powers of the mo- 
mentum p [49] , thus resulting in sextic and higher powers 
of p inside the square brackets of the effective Hamilto- 
nian (2.3). However, as shown in Ref. [49] and discussed 
in the Appendix, the mapping between the real and the 
effective descriptions can always be arranged in such a 
way that these "non-geodesic" terms are proportional to 
the radial momentum p r = n ■ p, thus vanishing for cir- 
cular orbits, at all PN orders. Moreover, we will show in 
the Appendix that this holds true even if the mapping 
(2.4) between the real and the effective Hamiltonians is 
assumed to be valid at all PN orders. This standpoint 
was also adopted in Ref. [79] , and it is the one that we em- 
brace in this paper. In other words, our results are valid 
for the class of EOB models that adopt the mapping (2.4) 
at all PN orders, and for which the non-geodesic higher- 
order momentum terms are proportional to p r , thus van- 
ishing for circular orbits. 

We notice that the 3PN-accurate expression (2.2a) of 
the EOB potential A{r) does not contain terms propor- 
tional to v 2 and v 3 , despite the fact the PN binding en- 
ergy does contain such terms. By contrast, at 3PN order, 
the EOB "radial" potential B{r) — and hence the inverse 
product D = (AB)^ 1 [cf. Eq. (2.2b)] — contains a term 
proportional to v 2 . The GSF results will allow us to con- 
trol the exact contributions proportional to the binary's 
mass ratio q = toi/to2 = v + 0{v 2 ) 3 thus only the terms 
linear in v in A(r) and D(r). 

In the next two Subsections, we will derive the GSF 
contributions to the EOB potentials, namely the func- 
tions A SF (u) and D SF (u) such that 4 

A{u) = 1-2u + isA S f(u)+0(v 2 ), (2.5a) 

D{u) = 1 + vD SF {u) +0{v 2 ). (2.5b) 

A. Self-force contribution to the EOB effective 
metric potential A 

We will restrict the discussion to circular orbits, com- 
puting first the EOB energy -Eeob for such orbits. Now, 



3 Without any loss of generality, we assume m\ ^ 7712- 

4 In the notations of Refs. [79, 81], we have Agp = a and Sgp = d. 
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the angular momentum L = can be determined as a 
function of the inverse separation u by solving the equa- 
tion 

Pr = ~ dH Q° B ( r ; Pr = Q,P4>)=0, (2.6) 

which is valid for circular orbits only. From the expres- 
sions (2.3) and (2.4) of the effective and EOB Hamilto- 
nians, this gives 



L 2 (u) = 

fj, 2 M 2 2uA(u) +u 2 A'(u) ' 



A'(u) 



(2.7) 



where we denote A' = dA/du. Replacing this result 
in the expressions for the Hamiltonians, we obtain the 
circular-orbit EOB energy as 



(2.8a) 



E eS (u) = (i. 



2A 2 (u) 



2A{u) +uA'{u) 



We also introduce the dimensionless coordinate-invariant 
quantity x = (MSI) 2 / 3 , where Q is the constant circular- 
orbit frequency, given by 



(2.9) 



as well as tq = M/x, a convenient invariant measure of 
the orbital separation. 

Since we are interested in improving the EOB Hamil- 
tonian using GSF results, we will work at linear order in 
the symmetric mass ratio v, thus neglecting terms 0(v 2 ) 
or higher. Inverting Eq. (2.9) and making use of the ex- 
pression (2.7) for the angular momentum yields [79] 



x ., , , 2 / l-2x 
6 SFV ' 3 



1 



Oiy 2 



(2.10) 

Next, we compute the (specific) EOB binding energy 
£"eob = (-^eob — M)//j> at linear order in v, for cir- 
cular orbits. From Eqs. (2.8), in which we substitute the 
coordinate-dependent variable u in favor of the invariant 
variable x using (2.10), we obtain 



■Eeob(x) 



1 - 2x 



1 



I -Ax A SF (x) 
(1 -3a;) 3 / 2 2 



x A' SF (x) 



vT 



3a; 
- 6x 



l-2x 
3 V VI -3a; 
If 1 - 2a; 



3(l-3x) 3 / 2 2\^T^3x~ 



- 1 



+ 0(u 2 ). (2.11) 



Recently, Ref. [87] used the first law of binary point- 
particle mechanics [88], together with GSF results for 
the redshift observable z [38, 82-85], to compute the 
binding energy E and total angular momentum L of 



a circular-orbit non-spinning compact binary, through 
next-to-leading order in the symmetric mass ratio v (at 
all PN orders). In particular, the specific binding energy 
E = E/n reads [87] 



E(x) = 



l-2x 
Vl - 3a; 



- 1 



z S f(x) 



Vl - 3.t - 1 



2 

x 7 - 24a; 
6 (1 - 3a;) 3 / 2 



-SF 



[x) 



+ 0(v 2 ), (2.12) 



where zsf(x) is the self- force contribution to the redshift 
z of the lightest point mass, which is known numerically, 
with high accuracy, for circular orbits down to rn = 5M 
(see below for more details). 

By construction of the EOB model, the EOB binding 
energy £"eob must coincide with the binding energy E 
of the real binary system when expanded in powers of v. 
Equating Eqs. (2.11) and (2.12) yields the following linear 
first-order ordinary differential equation for Asf(x): 



( 2 ' 8b ) 2xA' SF (x) 



1-Ax 
'l - 3a; 



A SF (a;) 



2a; Zg F (a;) — 3zsF(a;) + x ■ 



1 — 6x 

= x 

1 - 3a; 

5a; + 12a; 2 



(l-3a;) 5 



Vl - 3a; x 
(2.13) 



Interestingly, this differential equation can be solved an- 
alytically in terms of zgp(x). The solution is particularly 
simple, and explicitly reads 5 

A SF (x) = VT^ z SF (x)^xi [ l + - / ==j . (2.14) 

This is one of the most important results of this paper: we 
have succeeded in relating the known GSF contribution 
zsf(x) in the redshift observable to the EOB dynamics 
for circular orbits. As a result, we can now compute the 
EOB radial potential A(r) given in Eq. (2.5a) including 
all PN corrections, at linear order in v. As a consistency 
check, the PN expansion of 0sf(x), as given by Eq. (4.16) 
and Table I of Ref. [88], can be used with Eq. (2.14) to 
recover the 3PN expansion (2.2a) of Asf(x). 

We stress that the EOB model with this form of 
the potential A(r) automatically reproduces the shift of 
the ISCO frequency under the effect of the conserva- 
tive part of the GSF, as calculated by Barack and Sago 
in Refs. [39, 80]. This is because the notion of ISCO 
coincides with that of minimum energy circular orbit 
(MECO) for any orbital dynamics that can be derived 
from a Hamiltonian [89], and the MECO as computed 
from Eq. (2.12) returns the correct ISCO frequency shift 
[87]. This is a considerable improvement over earlier ver- 
sions of the EOB potential A(r), which were resorting 
to free parameters regulating unknown high-order PN 



5 The homogeneous solution must vanish because the PN expan- 
sion of Agp(x) cannot involve half-integer powers of x. 
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effects [53, 71, 79] in order to reproduce the fSCO fre- 
quency shift due to the conservative GSF. 

Furthermore, if we use our newly derived potential 
A(r) together with Eqs. (2.7) and (2.10) to compute the 
dimensionless angular momentum L = L/(pM) for circu- 
lar orbits, expressed in terms of the coordinate-invariant 
quantity x, we recover the result obtained in Ref. [87], 
namely 



L(x) 



3x) 
1 



1 



15x 



6 SF\ 



6^x~ (1 - 3x) 3 / 2 



+ 0(v 2 ), (2.15) 



which holds at all PN orders and at linear order in v. 
This comes at no surprise because the binding energy and 
total angular momentum for circular orbits satisfy the 
exact relation dE/dL — Mfl in the EOB model. There- 
fore, once the energy coincides with the exact expression 
established in Ref. [87], so must the angular momentum. 

The GSF contribution zsf(£) to the redshift observable 
has been calculated numerically in Refs. [38, 82-85] for 
a variety of orbital separations, in the range 5M ^ rn ^ 
500Af. In the companion paper [87], we established that 
this numerical data can conveniently be represented with 
an accuracy better than a part in 10 5 , using the compact 
analytical formula 



zsf(x) = 2x 



1 + CllX + CL2X 

1 + CL3X + a^x 2 + a$x 3 



(2.16) 



with the coefficients a x = -2.18522, a 2 = 1.05185, 
a 3 = -2.43395, a 4 = 0.400665, and a 5 = -5.9991. This 
fitting formula accounts for the leading-order (1PN) be- 
havior z SF (x) = 2x + 0(x 2 ) when x -> [88]. Since 
current GSF data for zc, F (x) is limited to orbital separa- 
tions ro ^ 5M, the GSF-accurate EOB potential A(u) 
given by Eqs. (2.5a) and (2.14) is (for now) only known 
in the range ^ x < 1/5. 



B. Self-force contribution to the EOB effective 
metric potential D 

The non-circular conservative dynamics of spinless bi- 
naries is regulated by the radial frequency Q r and by the 
averaged angular frequency fi^, respectively defined by 



Sir 



2tt 
P~ 
1 

P 



4>(t) dt = Kti r , 



(2.17a) 
(2.17b) 



where P is the radial period, namely the time interval 
between two successive periastron passages, 4> = dcj)/dt 
is the time derivative of the orbital phase </>(£), and 
A$/(27r) = K — 1 is the fractional advance of the pe- 
riastron per radial period. In the circular-orbit limit, 



by definition the radial frequency vanishes at the ISCO; 
hence the periastron advance K = n^/Qr blows up there. 
For this reason, Refs. [79, 81] found it convenient to 
work with the quantity W = 1/K 2 , which is regular at 
the ISCO. Reference [81] calculated numerically the GSF 
contribution to W, i.e. the function psf{x) such that 



W(x) = 1 - 6x + vp SF (x) + 0{v 2 ) . 



(2.18) 



The authors performed several fits of the GSF data for 
Psf{x) (in the range 6M ^ ^ 80M). In particular, 
they found that this data can be accurately reproduced 
at the 2.4 x 10~ 3 level by means of the compact analytic 
formula 



1 + 61 x 



I + x + b 3 x 2 



(2.19) 



with 61 = 13.3687, b 2 = 4.60958, and 63 = -9.47696. 
(Using a denser data set in a more limited frequency 
range, Ref. [90] later found that the values b\ — 12.9906, 
b 2 = 4.57724, and b 3 = -10.3124 yield a fit accurate at 
the 10" 5 level.) As with Eq. (2.16), the fitting formula 
(2.19) accounts for the leading-order (2PN) behavior of 
Psf(x) when x — > [81]. 

Reference [79] recently studied the dynamics of slightly 
eccentric orbits within the EOB formalism, and found 
that the GSF correction psf(^) to the periastron advance 
is related to the EOB potentials Asf(x) and Dsf(x) by 6 

/ 1 _ 2x \ 
PSF (x)=4xi [ l~ ^==J + A SF (x) + x A' SF (x) 

+ I (1 - 2x) A% F (x) + (1 - 6x) D SF (x) . (2.20) 

Solving the above equation for the unknown Dg F (x), we 
immediately obtain 



1 



1 — Qx 



y t ( l-2 x 
p SF(a;)+ 4^ 7 _-l 



A SF (x) -xA' SF (x)--(l- 2x)A'£ F (x) 



, (2.21) 



where j4sf(£) is given explicitly in terms of z$ F (x) by 
Eq. (2.14) above. Equation (2.21) is another important 
result of this paper: the EOB potential D(r) governing 
the non-radial motion, as given by Eq. (2.5b), is now 
known exactly at linear order in v 1 through the known 
GSF contributions zsf(^) and psf{x) to the redshift ob- 
servable and periastron advance. 

Note that the apparent pole at the Schwarzschild ISCO 
(x = 1/6) in Eq. (2.21) must be canceled out by a factor 
(1 — 6x) in the numerator, because the potential D{r) is 



6 This result requires that the quartic power in the radial mo- 
mentum p r inside the square brackets of Eq. (2.3) be neglected. 
This is correct in the limit of small eccentricity e<l, for which 
p r ~ e; see Ref. [79] for more details. 
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perfectly regular at the ISCO. This can be verified using 
the fit (2.16) for z$f(x), which is accurate to within 10~ D , 
together with the fit (2.19) for pg F (x), which is accurate 
to within 2.4 x 1CP 3 using the coefficients b±, 62 and 63 
from Ref. [81], and to within 1CP 5 , although in a more 
limited frequency range, using the coefficients 61, 62 and 
63 from Ref. [90]. Using these fits, the behavior of Dsp(x) 
near the ISCO is of the form D SF (x) = e/(x-l/6)+0(l), 
where the dimensionless parameter e is about 4 x 10~ 5 
with the coefficients 61, 62 and 63 from Ref. [81], and 
3.4 x 10~ 5 with the coefficients 61, b 2 and 63 from Ref. [90]. 
One can therefore argue that the value of e is comparable 
to the accuracy of the fits and thus compatible with zero. 

Finally, when using the known PN expansions of the 
GSF contributions psf{%) and zsf(x) to the redshift ob- 
servable and periastron advance (see e.g. Eqs. (5.32) 
and (5.39) of Ref. [79], and Eq. (4.16) and Table I of 



Ref. [88]), we recover, as expected, the 3PN expansion 
(2.2b) of D SF (x). 



III. HIGH-ORDER POST-NEWTONIAN TERMS 
IN THE EOB POTENTIALS 

Recently, the authors of Ref. [88] derived a first law of 
mechanics for non-spinning compact objects modelled as 
point particles, and moving along exactly circular orbits. 
By making use of previous GSF results for the redshift 
observable [38, 84, 85], they could determine the numer- 
ical values of some previously unknown coefficients en- 
tering the expression of the circular-orbit binding energy 
E. Through 6PN order, the (specific) binding energy 
explicitly reads 



E(x) 
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27 
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35 
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3969 
128 
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24 





( 675 


)x 2 + 
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34445 205 



576 



35 



15 



v In x ) x + 



v 64(f) - 

264627 
1024 



448 
l5" 



v lnx ) x 



96 
45927 
512 



+ ve e (v) + v ef(v) In a; ) x h + o{x°) 



(3.1) 



r 



where the numerical values of the 4PN, 5PN and 6PN 
coefficients e 4 (0), e 5 (0), e 6 (0), and eg n (0) were found to 
be [88] 



e 4 (0) = +153.8803(1) , (3.2a) 

e 5 (0) = -55.13(3), (3.2b) 

e 6 (0) = +588(7), (3.2c) 

4 n (0) = -1144(2). (3.2d) 



(The uncertainty in the last digit is indicated in paren- 
thesis.) Note that the leading-order 4PN and next-to- 
leading order 5PN logarithmic contributions to the bind- 
ing energy are known analytically [85] . The value of the 
post-GSF coefficient (term cx v 2 ) in the 5PN logarithmic 
contribution is that corresponding to the "physical prob- 
lem" in the language of Ref. [88], i.e. when the helical 
symmetry is not imposed (see [88] for more details). 

In the next two Subsections, we will use Eqs. (3.1) and 
(3.2), together with the results of Ref. [81], to compute 
the 4PN, 5PN and 6PN coefficients in A(u), as well as 
the 4PN and 5PN coefficients in D(u), at linear order in 
the symmetric mass ratio v. 



A. High-order post-Newtonian terms in the EOB 
metric potential A 

Similarly to what was done in Sec. II, we restrict here 
to the class of EOB models that adopt the mapping (2.4) 
at all PN orders, and for which the non-geodesic higher- 
order momentum terms are proportional to the radial 
momentum. As a consequence, the 4PN, 5PN and 6PN 
corrections determined in Eqs. (3.1) and (3.2) enter the 
effective Hamiltonian only through the radial potential 
A(r), i.e., through the coefficients a§{v), o}§(y), oq{v), 
ajf(f), a-j(v) and a}f(v ) in the PN expansion 

, , . „ 3 /94 41tt 2 \ 4 

A(u) = 1 - 2u + 2vu 3 + ( — - — i/w 4 

+ v [a 5 (v) + 4 n (f) law] u 5 
+ u [a 6 (f) + 4 a (z/) law] u 6 
+ v [a 7 {v) + 4 n (i/) Inw] u 7 + o(u 7 ) . (3.3) 

A comment regarding the general structure of the PN 
expansion of A(u) is in order. It was argued in Ref. [85] 
that terms involving powers of logarithms should not oc- 
cur in the conservative part of the dynamics of a compact 
binary system before the very high 7PN order. Within 
the class of EOB Hamiltonians that we consider in this 
work, the potential A is directly related to the conserved 
binding energy E of the real binary system. Hence we 



7 



expect no term of the type (lnu) p , with p ^ 2, in the 
6PN-accurate expansion (3.3) of A(u). 

Now, in order to derive the v — » limits of the coef- 
ficients ai{v) and a' n (f) (with i = 5,6,7), we first de- 
termine the angular momentum L — p^ as a function of 
the circular-orbit frequency fi by solving Eqs. (2.6) and 
(2.9). Second, we insert L = L(fi) in Eqs. (2.3) and (2.4), 
set p r = 0, and expand Eq. (2.4) in PN orders. Third, 
after subtracting the rest-mass contribution, we equate 
the result to Eq. (3.1), and finally obtain 
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C(o). 



(3.4a) 
(3.4b) 

(3.4c) 
(3.4d) 

(3.4e) 
(3.4f) 



Note that the values of the 4PN and 5PN log-coefficients 
ag n (0) and a{f(0) are known analytically. 7 Using the nu- 
merical results (3.2) for the binding energy, we find for 
the other 4PN, 5PN and 6PN unknown coefficients: 



a 5 (0) = +23.50190(5) , 
a 6 (Q) = -131.72(1), 
a 7 (0) = +118(2), 
4 n (0) = -255.0(5) . 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 



As a consistency check, we verified that the same values 
for these coefficients also follow from expanding in PN 
orders the potential A given by Eqs. (2.5a) and (2.14), 
using the 6PN-accurate fit to zsf(x) given in Eq. (4.16) 
and Table II of Ref. [88]. 

Finally, we note that the values of the coefficients a 5 (0) 
and ae(0), as determined in Ref. [79] [see Eq. (4.46) there] 
by combining the EOB /NR comparison of Ref. [70] with 
the constraint coming from the GSF calculation of the 
ISCO frequency shift, are in poor agreement with the 
exact results (3.5a) and (3.5b); in particular the signs 
of both coefficients are wrongly predicted. Even with- 
out taking into account the constraint from the GSF, 
the correct values (3.5a) and (3.5b) lie far outside the 
(banana-shaped) region of the 05(0), ae(0) plane favored 



7 Actually these logarithmic terms are known for all mass ratios: 
4» = f and «*(„) = _ lf u . 



by a comparison of the EOB prediction with the NR 
waveform of an equal-mass binary black-hole simulation 
[70, 91]. This disagreement does not surprise us. Indeed, 
several causes can affect the correct values of the coef- 
ficients 05(0) and ag(0) when they are extracted from a 
EOB/NR calibration, such as radiation-reaction effects 
and all higher-order PN contributions in A(r), which be- 
come important during the plunge and close to merger, 
but are neglected in those calibrations. 



B. High-order post-Newtonian terms in the EOB 
metric potential D 

We will now show that by combining the results of 
this paper with those of Ref. [81], we can also deter- 
mine the exact or approximate numerical values of the 
non- logarithmic and logarithmic coefficients c/4 (z/), ds(y), 
d^{v) and d^(v) appearing at 4PN and 5PN order in the 
potential D: 

D(u) = 1 + 6v u 2 + (52 - 6v) v u 3 

+ v[d i {v) + d}£{v) lnw] m 4 

+ v \d 5 (u) + 4 n H lnu ] " 5 + °( u5 ) ■ 



(3.6) 



The authors of Ref. [81] extracted numerically some un- 
known high-order PN terms in the function pgp(x), and 
used the relation (2.20) between psf, A§f, and D$f to 
put constraints on some unknown high-order PN coeffi- 
cients appearing in the EOB potentials A(u) and D(u). 
In particular, Ref. [81] derived the following constraints 
on the v — > limit of the PN coefficients d±(v) and d§(v): 



10o s (0) + d 4 (0) + -4 n (°) = 518.6tJ, (3.7a) 



14o 5 (0) + 6d 4 (0) - 15a 6 (0) -d 5 (0) 
+ 84 n (0)-ya 6 n (0)= 4779;^ , 



(3.7b) 



as well as the following exact relations between the v — » 
limit of the PN coefficients d^{v) and d}£(v): 



10 «#(()) + <*5[ , (o) 



2 512 
15 



14 4 n (0) - 15a 6 n (0) +6< (0) - 4"(0) 



(3.8a) 
11336 
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(3.8b) 



Substituting Eqs. (3.4b), (3.4d), and (3.5) in the above 
equations, we get 



d 4 (0) = +226.0+; , 
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(3.9a) 
(3.9b) 
(3.9c) 
(3.9d) 
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FIG. 1: The fractional difference between the GSF-accurate (specific) binding energy and the (specific) binding energy calculated 
at linear order in v using standard PN approximants (upper plot), and EOB approximants (lower plot), for v — 0.01 (left panel) 
and v — 0.1 (right panel). 
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FIG. 2: Same as in Fig. 1, for an equal-mass binary (y = 0.25). 



As a consistency check, we verified that the values (3.9) 
can be recovered by inserting the PN fits to Zsf(#) and 
Psf(^) given in Refs. [81, 88] into the exact expression 
(2.5b) and (2.21) for D. 

In summary, building on the results of Refs. [81, 88], 
we have computed the 4PN, 5PN and 6PN terms in the 
EOB potential A(u), as well as the 4PN and 5PN terms 
in the potential D(u), at linear order in v. 



IV. DISCUSSION AND CONCLUSIONS 

The potential A(u) given by Eqs. (2.5a) and (2.14) en- 
sures that the EOB binding energy for circular orbits co- 
incides, at linear order in v, with the exact result (2.12). 
In order to investigate the properties of the EOB resum- 
mation, as opposed to the standard PN expansion, we 
will here compare this exact result for the binding energy 
to the PN predictions, as well as to the EOB binding en- 
ergy, as computed with the PN-expanded version of the 
potential A(u) [cf. Eq. (3.3)]. 

In Figs. 1 and 2 we show the fractional difference be- 
tween the GSF-accurate (specific) binding energy E(x) 
given by Eq. (2.12), and either the EOB (specific) binding 
energy obtained from Eqs. (2.8) using the PN-expanded 
potential A (at linear order in v), or the PN (specific) 
binding energy given in Eq. (3.1), including only the 
terms linear in v. The fractional difference is presented 
as a function of x = (A/SI) 2 / 3 up to x = 1/5. We con- 
sider three representative symmetric mass ratios, namely 
v = 0.01,0.1, and 0.25, and several PN orders. 8 

From Figs. 1 and 2 we observe that up to a; = 1/5 
the EOB-approximants are much closer than the PN- 
approximants to the exact GSF result for the small mass- 
ratio case v — 0.01, and are (roughly) comparable, with 
some differences depending on the PN order, to the exact 



8 In order to express the EOB binding energy (2.8) at nPN order 
as a function of x rather than u, we insert Eq. (2.7) into the 
relation £ 3 / 2 = MdHEOB/dL(r = l/u,p r = 0, L), and invert 
the latter to obtain u as a PN expansion in x. 
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GSF result for the comparable mass-ratio cases v = 0.1 
and v — 0.25. 

These results confirm the utility of resumming the PN 
dynamics around the test-mass limit if the goal is to ob- 
tain a flexible model that can smoothly bridge between 
the test-mass and equal-mass limits. This flexibility ex- 
ists in the EOB model not only at the level of the bind- 
ing energy for circular orbits, but more importantly at 
the level of the Hamiltonian, thus for generic orbits and 
beyond the innermost stable circular orbit. It is a cru- 
cial feature that has allowed to build faithful inspiral- 
merger-ringdown templates that can span the entire bi- 
nary's mass-ratio range [57, 60-78]. 

We leave to future work a detailed study of how the 
EOB model, augmented with the gravitational self-force 
results in A(u) and D(u), performs against numerical- 
relativity simulations of comparable-mass black-hole bi- 
naries [71]. Comparisons using a 4PN, 5PN, or even 6PN- 
accurate EOB model can already be pursued. However, 
in order to use the EOB potentials with all PN terms 
linear in v, more self-force data is needed for the redshift 
observable z§y(x) beyond x — 1/5, and for the periastron 
advance psf(x) beyond x = 1/6. 
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Appendix A: On the generic structure of the 
effective-one-body Hamiltonian 

Following Ref. [49], but including also the logarithmic 
contributions that were left out of their analysis, we show 
that it is possible to build EOB Hamiltonians such that 
(i) the mapping (2.4) holds at all PN orders and (ii) the 
non-geodesic terms in the square brackets of Eq. (2.3) are 
proportional to the radial momentum, at all PN orders, 
thus vanishing for circular orbits. 

We start by considering the generic structure of the 
(specific) two-body Hamiltonian in the center-of-mass 
frame, at a given PN order, as a function of the reduced 
canonical variables p = p/p, and q = r/M : 



H (n , k)PN (r,p) = (lnq) k |p 2 (" +1 ) 
1 



+ - \p 2n +p 2n - 2 (np) 2 + ■■■ + (np) 2n ] 
Q 
1 



2(n-l) 



(np) 



+ ■ 



(Al) 



where we introduce the notations q = ^/q ■ q and (np) = 
n ■ p, and where we use a subscript (n, fc)PN to denote 
the contribution to the nPN-accurate Hamiltonian which 
is 0(c~ 2n ) and proportional to (ln<j) fe , with k ^ 0. In- 
deed, the general structure of the near-zone expansion 
(formally r/c — > 0) of the PN metric is known to be 
of the type [In (r / c)] k (r / c) n [92], yielding terms propor- 
tional to (lnq) k in the conservative dynamics starting at 
4PN order [93] (see also [85] for a more recent discussion) . 

The number of independent coefficients in the Hamil- 
tonian (Al) is easily found to be [49] 



(n + 1) (n + 2) 



1 



(A2) 



The mapping of the real Hamiltonian onto the effective 
Hamiltonian can be achieved through a canonical trans- 
formation [47, 49] . The most general generating function 
reads 



G, 



ni)PN 



(r,p) = (\nq) k (r ■ p) lp 



~2n 



l r 

(1 



2(n-l) 



(np) 



2(n-l) 



(A3) 



Thus, the number of independent coefficients in the gen- 
erating function is [49] 



C G (n,k) 



n (n + 1) 



1 . 



(A4) 



At 3PN or higher orders, non-geodesic (NG) terms result- 
ing in quartic or higher powers of the momentum p can 
appear inside the square brackets of the effective Hamil- 
tonian (2.3). At 3PN order, the NG terms symbolically 
read 



Q 3 p N (r,p) = ^ [P 4 +f ( n P) 2 
q A L 



(np) 4 



(A5) 



(The first logarithms appear at 4PN order.) At higher 
PN orders (n > 3), the generic structure is 



Q(n,k)PN(r,p) = (hig) fc 




(np) 2 



^ [p 4 +f (np) 2 + (np) 4 



(A6) 



It is straightforward to derive that the number of arbi- 
trary coefficients in Q( n ,fc)pN is 



C NG (n,k) 



(n + 3) (n - 2) 



(A7) 



hence the number of NG terms that depend on the radial 
momentum (np) is 



C^(n,k) = C NG (n,k)-(n-2) 



(n + l)(n-2) 



(A8) 
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Given that at the PN order (n, k) we have two new coef- 
ficients in the effective metric potentials A and D, mul- 
tiplied respectively by u n+1 (\nu) k and u n (\nu) k , but no 
new coefficient in the mapping between the real and effec- 
tive Hamiltonians, because we assume that (2.4) is valid 
at all PN orders, we obtain that the difference between 
the number of equations to satisfy and the number of 
unknowns is 

S(n, k) = C H (n, k) - C G (n, k) - 2 - C^(n, k) 



which is always zero or negative starting at 3PN (n = 3). 
This proves that for all n ^ 3 (and any k J? 0), it is possi- 
ble to build the mapping between the real and the effec- 
tive descriptions in such a way that Eq. (2.4) holds, and 
circular orbits follow from a "geodesic" effective Hamil- 
tonian [i.e., all the higher-momentum terms inside the 
square brackets in Eq. (2.3) can be chosen to be propor- 
tional to the radial momentum (np)]. It is for this class 
of EOB Hamiltonians that we have determined the effec- 
tive metric potentials A and D at all PN orders, linearly 
in the symmetric mass ratio v. 
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